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ABSTRACT 

HD 179821 is an evolved star of unknown progenitor mass range (either post- 
Asymptotic Giant Branch or post-Red Supergiant) exhibiting a double peaked spectral 
energy distribution (SED) with a sharp rise from ~ 8 — 20 /xm. Such features have been 
associated with ejected dust shells or inwardly truncated circumstellar discs. In order to 
compare SEDs from both systems, we employ a spherically symmetric radiative transfer 
code and compare it to a radiative, inwardly truncated disc code. As a case study, we 
model the broad-band SED of HD 179821 using both codes. Shortward of 40 ;um, we 
find that both models produce equivalent fits to the data. However, longward of 40 /im, 
the radial density distribution and corresponding broad range of disc temperatures pro- 
duce excess emission above our spherically symmetric solutions and the observations. 
For HD 179821, our best fit consists of a T^ff = 7000 K central source characterized by 
TV ~ 1-95 and surrounded by a radiatively driven, spherically symmetric dust shell. The 
extinction of the central source reddens the broad-band colours so that they resemble 
a Teff = 5750 K photosphere. We believe that HD 179821 contains a hotter central 
star than previously thought. Our results provide an initial step towards a technique 
to distinguish geometric differences from spectral modeling. 

Subject headings: stars: AGB and post-AGB - radiative transfer - planetary nebulae: 
general - infrared: general - circumstellar matter - submillimetre 



- 2 - 



Introduction 



HD 179821 (IRAS 19114+0002; V1427 Aql; SAO 124414; AFGL 2423) is an evolved star 
surrounded by gas and dust ejected during a phase of mass-loss. The lumi nosity of this object is 
undetermined as the Hipparcos parallax measurement ([0.18 it 1.12 mas]; iPerrvman et al.l 119971 ) 
allows for any distance greater than 1 kpc. HD 179821 is either a close post-Asymptotic Giant 
Branch (post-AGB) star {D = 1 kpc, Mi ~ 3 — 4 Mq) or a distant, massive (-D = 5 — 6 kpc, 
Mj ~ 20 — 30 Mq) post-Red Supergiant (post-RSG). The post-AGB phase of stellar evolution is 
a short ('^ 10^ yrs) period in which initially intermediate-mass main sequence stars (< 8 Mq) 
transition from the Asymptotic Giant Branch to the planetary nebula (PN) phase. In contrast, 
the progenitors of Red Supergiants are massive main-sequence stars (> 8 Mq) evolving toward the 
Wolf-Rayet stage and may end their l i ves as supernovae. For a short review on the post-RSG status 
of HD 179821, see lOudmaiier et al. I (|2008l ). In both phases, mass is ejected into the circumstellar 



environment. Distinguishing between these classes of objects can be difficult. 

In this paper, we model the broad-band SED of HD 179821 using two distinct radiative transfer 
codes. The first code, "DUSTY", computes rad iative transfer th rough a spherically symmetric 



shell with canonical density profile p{r) cx r ^ ( Ivezic et al.lll999l ). The second code computes 
radiative transfer through a n axisymmetric, flared disc with canonical density profile p{r) oc r~^/^ 



( D'Alessio et al. I I2OO1I .I I2OO5I ). While we believe the mid-IR and radio imaging of HD 179821 indicates 
a roughly spherical nebula, the point of this paper is to determine whether spectral modeling, 
through the use of distinct radiative transfer codes, can constrain geometric features of evolved 
star nebulae. Thus, HD 179821 serves as a test case for this investigation. 

We model HD 179821 for both nebular geometries and investigate the degeneracies between 
the SED model fits. In addition, we present and compare the ISO spectrum of HD 179821 to the 
post-AGB object HD 161796 and post-RSG object IRC +10420. All three objects display a strong 
mid-IR excess and exhibit steep increases in their SED's beyond 8 microns. In Section 2, we discuss 
the enigmatic nature of HD 179821 in both the post-RSG and post-AGB scenarios. In Section 3, 
we describe our synthetic photospheres and extinction corrections. In Section 4, we present our 
SED model fits from the radiative transfer disc code. In Section 5, we present our model fits from 
DUSTY and identify wavelength regions in which degeneracies occur between the two codes. 

For intermediate mass stars below 8 Mq, the late stages of stellar evolution are characterized 
by transitions from quasi-spherical mass- losing AGB stars to complex aspherical post-AGB objects. 
Many post-AGB objects feature toroidal density enhancements, bipola r jets and an array of other 
non-spherically symmetric features, the origin of which is unknown (jBalick &: Frank! |2002| ). In 
addition, most post -AGB objects display a large momentum excess above what would be supplied by 
radiation pressure (jBujarrabal et al.ll200ll ). While the origin of the additional source of momentum 



is unclear, a binary companion is an attractiv e candidate as energy and an gular momentum can 



transfer from the secondary to th e primarv (INo rdhaus B 



supported by recent observational (jDe Marco et al.,, 2004 . 120071 ) and theoretical (jMoe &: De Marco 



ackmani l2006l ). This hypothesis is 
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20061 : ISokeij 120061 ) efforts suggesting that a significant fraction of PNe may be descendants of 
interacting binaries. Binary companions can influence mass loss in many ways. A common envelope 
phase, even for lo w-mass companions, can lead to equatorial outflows, poloidal outflows and the 
formation of discs (IReyes-Ruiz Lopea Il999l . iNordhaus &: BlackmanI 120061 . iNordhaus et al.l 120071 . 
Nordhaus &: BlackmanI 120071 ) . For wider binaries, low-mass c ompanions can iii duce spiral waves 
and convert amorphous dust to crystalline dust via annealing (jEdgar et al.l 120071 ). 



For massive stars, the end stages of stellar evolution are poorly understood. The post-red 
supergiant (post-RSG) phase is among the most luminous and uncertain epochs of post-main 
sequence stellar evolution. Like their interme diate mass conter parts, post-RSG's are thought to be 
ejecting a substantial portion of their mass (|Decin et al.l 120061 ). However, a massive circumstellar 
envelope has only been detected (scattered light, IR emission, molecular lin e emission) around the 
post- RSG IRC +10420, thus making a comparison with HD 179821 difficult (jKastner Weintraub 
199d ). 



Mass-loss may be consistent with a radiation driven spherically symmetric outflow (|Castro-Carrizo et aL 
20071). However, if H P 179821 does indeed display a momentum excess in its ejected nebula 



(jBujarrabal et al.ll200ll ). then it is not unreasonable for a binary companion to have influence d the 
outflow and produced asymmetries. Interestingly, the WFPC2 images of lUeta et al.l l2000l show 
collimated bipolar structures emerging from the dust shell. These may be remnant jets propagat- 
ing through the envelope. There is a l so evidence for slight clump y regions in both OH maser and 



mid-IR emission (|Gledhill et ahlboOll : Ix^stner fc Wein traubll l99d). A density asymmetry may also 



be present as the ^^CO line profiles are asymmetric (jBujarr abal et al.lll992l ). 



The spectral energy distribution of HD 179821 exhibits a double-peaked shape, indicative of 
a stellar photospheric component and ejected dusty component. The spectral energy distribution 
(SED) is consistent with photospheric emission to ~ 8 ^m at which point there is a steep increase 
until the second peak at ~ 25 /im. The overall SED of HD 179821 is t hus remarkably simila r 
to that of the transitional young stellar object CoKu Tau/4 (see Fig. 2 of lD'Alessio et al.ll2005l ). 
which has been well modeled as a dusty disc with an inner hole. The sharp, interior wall of the 
disc is illuminated by the central source and produces the excess infrared emission. HD 179821 
also exhibits an evacuated interior region between the central object and hence, the system might 
be modeled similarly to CoKu Tau/4. However, the dust envelope of HD 179821 could also be 



a detached spherical sheU and past modeling effort s have treated it as such (ISurendiranath et al 



2002 



Buemi et all 120071 : ICastro-Carrizo et al.l 120071 ). 



2. HD 179821: post-AGB or red supergiant? 



HD 179821 (ga lactic coordinates (/ = 35.62°, b =-4.96°) is usually classified as a G5 star 
(jHrivnak et al.lll989l : see §2.1). It exhibits large mm- w ave CO line widths (^ 70 km/s) centered on 
the local standard of rest velocity Vlsr ~ 100 km/s (jZuckerman &: Dyckl 119861 : iBujarrabal et al. 
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19921 ). OH masers have been detected (jLikkellll989l : iGledhill et alj|200ll ) and the presence of a ^ 10 
sihcate feature confirm that this object is oxygen-rich (0/C ~ 2.6: lReddv &: Hrivnaklll999l ). 



On the whole, the chemical composition of H P 179821 differs from th at of an average F 
supergiant and is consistent with a post-AGB star (jReddy &: Hrivnald Il999l ). In particular, an 
overabundance of s-process elements suggest extra-mixing during the AGB phase and is further 
supported by a low isotopic ^^C/^^C < 5 ratio consistent with deep mixing in a post-AGB object 



(jCharbonnellll995l : iJosselin &: Lebrdl200lh . Underabundances of carbon and the s-process element 
zirconium suggest that HD 179821 left the AGB before the third dredge up, consistent with its 
0-rich properties. The presence of an active photochemistry and the p roduction of HCO"*" ar e 



detected in the outflows of many post-AGB objects including HD 179821 (IJosselin &: Lebre 



Currently HCO"'" has only been detected in one red supergiant star (VY CMa; IZiurys et al 



2001 ) 



20071 ). 



In general, the chemical composition is similar to that of high latitude , hot post-AGB stars and 



in pa rticular bears a strong resemblance to HD 161796 (jThevenin et al.l l2000l : iReddv Sz Hrivnak 



19991 ). The ISO SWS spectrum of HD 179821 and broad-band SED also appear very similar to that 
of the post-AGB star HD161796 (= IRAS 17436-^5003 = V814 Her = SAO 30548) (see Fig. d]). 
Also shown in Fig. fflis the ISO SWS spectrum of the post-RSG IRC -M0420 (=V1302 Aql = IRAS 
19244-1-1115). All three spectra show a substantial infrared excess, however there are differences 
between HD 161796 and IRC -^10420. The similarities in the ISO spectra for HD 179821 and HD 
161796 suggest that both systems have undergone similar mass-loss histories. 

Based on the kinematic distance of ~ 6 kpc, HD 179821 is a post-RSG. However, this result 
should be viewed with caution as the meth ods used to derive the kinematic distance are not valid 
for the galactic coordinates of HD 179821 (jjossehn k LebrelboOlh . If in fact HD 179821 is 6 kpc 
away, then its latitude indicates a position of 500 pc a bove the galactic plane. T his location is 
approximately 5-6 times the scale height for supergiants (IScheffler &: Elsaesseiill987l ). 



The strength of the 01 triplet line and the distance deter mined from the interste llar Na I 
Dl and D2 components imply that the object is very luminous (jReddv &: Hrivnaklll999l). Na and 
Na I emission at 2.21 fim is similar to that detected in known supergiants (jHrivnak et al.lll994l ) 
and abundances of sor t ie mo lecules are 40 times lower than the averages in post-AGB objects 
( Quintana-Lacaci et al.| 2007). CO outflow observations also support the post-RSG hypothesis 



Bujarrabal et al 



2001 



). In particular, they found a quasi-spherical, unique component to the 
wind with a large outflow velocity of Vexp = 33 — 35 km/s. This is higher tha n most post-AGB 



spherical wind outflow velocities (~ 10 — 15 km/s) (iKastner fc Weintraublll995l ). 



In short, whether this object is a post-AGB star evolving toward a PNe phase or a post-RSG 
star whose fate is to explode as a SNe remains unresolved. 
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3. Photospheric Models and Extinction 



Because of uncertainty in line-of-siglit visual and IR extinction (interstellar + circumstellar), 
the reduction in magnitude of the HD 179821 spectrum at a given wavelength is unknown. The 
inferred extinction must at least be enough to make the maxima of t he double peaks equal on 
the de-reddened SEP. Adopted values of A; have ranged from 1.8 to 4 (Ivan der Veen et al J 1 1994 : 
Hrivnak et alJll989l : ISurendiranath et alJl2002l : iHawkins et al.lll995l ). 



The choice of can greatly influence the inferred geometry of the system. The flux ratio of the 
peak of the infrared excess to that of the peak of the photosphere, a, determines the fraction of the 
central radiation intercepted by circumstellar material. If a ~ 1, then the circumstellar material 
is intercepting most of the radiation from the central object, corresponding to a spherical shell 
geometry. If, however, a < 1, the geometry can be torus- or disc-like with radiation from the central 
source escaping through bipolar cavities. Increasing raises the corrected photospheric emission 
relative to the longer infrared wavelengths, which are only marginally affected by extinction. The 
photometric data is presented in Tables 1 and 2. 

Since the extinction is unknown, the temperature of the central source may be uncertain. 
Low- resolutio n spectra hav e classi fied HD 179821 as a l uminous G-supergia. nt (G la iBidelman 



1981 



G4 0-Ia lKeenanlll983l : G5 Ia lBuscombelll984l : G5 la lHrivnak et al.l 11989 ) while the spectral 



ty pe inferred from high-res olution spectra indicates T^ff = 6800 K ( Zacs et al. 199q'). T^f f = 6750 
K (jReddy &: Hrivnaklll999l ) corresponding to an F star. In addition, lArkhipova et al.ll200ll conclude 
from their 1999 spectra that HD 179821 was an F star in the 1990's. If the extinction is purely 
interstellar, then a synthetic photosphere of Tg// = 5750 K (G spectral type) would correspo nd to 



A^ = 2 while Tg// = 7000 K (F spectral type) would correspond to A^ = 3 (|Kurucall993l ). For 
both radiative transfer models, we assume the central source is a Tf^jf = 7000 K star. However 
if the surrounding nebula is optically thick, the broad-band colours are reddened and the central 
source obscured. 

We model the full spectral energy distribution of HD 179821 in order to fit photometric (optical, 
2MASS, sub-millimeter, millimeter) and spectroscopic (ISO) data. In particular, we compare SED 
model fits from both codes to investigate degeneracies between circumstellar dust geometries. In 
Section 4, we present an inwardly truncated, flared disc model. The inner wall absorbs radiation 
from the central source and marks the transition to the optically thick outer disc. The circumstellar 
material in this case corresponds to a torus- like geometry. In Section 5, we present a spherically 
symmetric shell model using the radiative transfer code DUSTY. 



4. Inner Wall, Edge-on Disc Models 



The central star of HD 179821 is assumed to be located at 5 kpc with the following properties: 
Te// = 7000 K, Ri, = 287 Rq, = 8 Mq. For these parameters, we are implicitly assuming that 
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the system is a post red supergiant. However, since our model scales with distance, our modeling 
results could also apply to a post-AGB system at 1 kpc. 

We consider an optically-thick disc with an inner wall at radius r = and height H^. The 
wall marks the transition between the evacuated interior region and the optically-thick, outer disc. 
The wall is illuminated at normal incidence by the central object and is assumed to be uni form 
in the vertical direction. The radial structure is solved according to iD'Alessio et al.l |2005| . At 



wavelengths > 35 //m, emission is dominated by the outer disc. Here, we focus on fitting the 
ISO spectrum in the case of an edge-on disc and discuss fitting the sub-millimeter and millimeter 
wavelengths in Section 5. 



3.5 



We assume spherical dust grains in the wall at mosphere of size d istribution n{a) ~ a 
between minimum and maximum radii amin and Umax (IMathis et al.lll977l ) . Opacities are calculated 



using Mie theory (jWiscombd Il979l ) . 



We consider several silicate compositions (jHenning et al 



1999 



Jager et al.l l2003l ) . For each 



composition, we vary two parameters not constrained a priori: height of the wall and disc inclination 



angle. In addition, we vary amin and amax 
index fixed. 



to obtain the best spectral fit, leaving the power-law 



The shape of the spectrum between ~ 8 — 30 /xm allows us to rule out several silicate composi- 
tions. The 10 /im silicate feature is smooth and does does not show evidence of the substructures 
associated with crystalline components. This most likely indicates that the silicates are glassy and 
amorphous. Therefore, we can im mediately disregar d crystalline silicates of mean cosmic composi- 
tion (Mgn cjFen A'liCan n'^Aln n^SiO'ii njaeger et al.lll994l ) as they contribute substructures not o bserved 
in the smooth silicate features. Additionally, the glassy silicates from I Jaeger et al.l Il994l . do not 
fit the spectrum longw ard of 30 nm and hence, we rule them out. Hot shell circumstellar silicates 
(jOssenkopf et al.lll992l ) do not provide a reasonab le fit to the ISO spectrum. If instead we use cool 
shell circumstellar silicates (jOssenkopf et al.lll992l ). we can fit the 10 /um feature but can not match 
the emission longward of 20 fj,m. 

We obtain reasonable fits to the ISO spe ctrum of HP 179821 u sing glassy pyroxene (optical 



Dorschner et a 



constants fromlHenning Sz Mutschkdll997l and 
Dorschner et al.l Il995l ) and glassy bronzite ( Dorschner et al 



19951), glassy olivine (MgiFeiSi04; 



19881 ). Fig. [2] presents our fits for 
glassy pyroxene for two different inclinations: fi = cos{i) = 0.25 (top) and /U = 0.4 5 (bottom). For 
both figures, the temperature, Tq, at the innermost radius of the wall is 128 K (see ID'Alessio et al 



20051 ). In addition, the minimum and maximum grain radii which best fit the ISO spectrum are 
given by amin = 0.005 fim and amax = 1-0 /im. We have tried many different values of a^ax over 
our range of parameters and grains of ~ 1 fim and larger are required for our best fits. For the top 
model in Fig. [21 the position of the wall is R^, = 3.7 x 10^ AU with a height of Hy, = 1.73 x 10^ AU. 
We find good agreement up to 20 fim at which point our model undershoots the observed SED. 

If fi is increased (i.e. tilted more toward face-on), as in the bottom figure, it is possible to 
increase the emission slightly at 20 fiui however the flux at 18 ^m starts to deviate from the ISO 
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spectrum. The position of the wall in this model is the same, but a slightly shorter wall is required; 
Hyj = 1.69 X 10^ AU. 

Fig. [3l compares disc models obtained for two different silicate compositions: glassy olivine 
(top) and glassy bronzite (bottom). For the top figure, the best fit requires a slightly hotter wall, 
T{ri) = 131 K, located at = 6.3 x 10^ AU. The waU height is given as H^j = 3.4 x 10^ AU. For 
olivine, we require amm = 0.005 //m and amax = 5.0 /xm for our best fit, however, while we find 
good overall agreement between 2 and 45 /im, our model produces excess emission between 15 and 
18 ^m. For bronzite, we find that Tq = 128 K with ii^ = 6.7 x 10^ AU and = 3.7 x 10^ AU. 
For this model, the grain sizes are amin = 0.005 ^m and amax = 2.0 /xm. 

Using all three silicate compositions, we find that = 5.2 it 1.5 x 10^ AU and Hyj = 
2.7 ± 1.0 X 10^ AU with the wall dust temperature ranging from 128 to 131 K. With a limited 
library of dust opacity data, we cannot determine the exact composition of the wall, other than 
to suggest that it could be a mixt ure of our suggested rn aterials or perhaps could be fit by other 



amorphous silicates. For instance. lJusttanont et al.lll992l fit the 10 and 20 ^m features from their 



UKIRT CGS3 spectra using a colder (~ 90 — 95 K) mixture of olivine and magnetite. 

To summarize, we found satisfactory fits to the ISO spectrum with the disc code but could not 
fit photometry longward of 40 /im. We aim to compare modeling results from spherical and disc 
codes in regimes which might produce overlap. Thus, we study the shell model in the next section. 



5. Spherical Shell Models 



We use the spherically symmetric, radiative transfer code DUSTY (jivezic et al.lll999l ). Given 
an incident radiation field and specified dusty region, the code self-consistenly calculates the emer- 
gent flux including dust scattering, absorption and emission in a spherical, non-rotating environ- 
ment. In addition, for radiatively driven winds, DUSTY computes wind structure by jointly solving 
the hydrodynamic equations. The dust density profile within the shell is given as p{r) oc r~^, con- 
sistent with a constant mass loss rate. 



DUSTY was previously used to model HD 179821. ISurendiranath et al.ll2002l produced a best 
fit shell model with an inner radius of Ti = 8.7 x 10^ AU and temperature T{ri) = 130 K. However, it 
is difficult to tell how good a fit it is as they omit the ISO spectra in their model. iBuemi et al.ll2007 
also modeled HD 179821 using DUSTY and obtained a slightly cooler shell {T{ri) = llOK) located 
at a farther distance (rj = 12.0 x 10^ AU). However, neither investigation used the ISO spectrum 
to constrain their models. We reproduced both previous DUSTY models. However, when overlaid 
on the SED with the added constraint of the ISO spectrum, neither model proved satisfactory. As 
the ISO spectrum provides the most stringent constraints on model parameters, our fit improves 
previous work. 



For our spherical model, the temperature at the inner shell boundary (r(rj) = 128 K) was 
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varied until the model infrared excess matched the ISO spectrum. We fix the grain distribution 



as n(a) oc a and vary the maximum, 



and minimum, a^mj grain sizes. The equilibrium 



temperature of a dust grain depends on size and composition. Small grains superheat while large 
grains exhibit blackbody emission. At each radius within the nebula, DUSTY assumes that the 
dust grains are at the same temperature. This approximation is a generic feature of both codes used 
in this paper and warrants acknowledgment. Additionally, both codes assume isotropic scattering. 
Future work should incorporate anisotropic scattering in both geometries. 

In addition, we vary the ratio of the outer to inner shell radius, Y = The central radiation 
source is assumed to be the same 7000 K model photosphere used previously in §4. The optical 
depth through the shell at 0.55 /im is also varied. We opt for a radiatively driven wind and vary 
the composition of the grains. 

Our best fit to the data yields an optical depth through the shell of Ty = 1-95 it 0.03. In 
addition, we find a good fit using minimum and maximum grain sizes amin = 0.005 //m and 
<imax = 0.25 yum. We find the ratio of the outer to inner shell radius produces similar results for 
F = ^ = 15±2. 



The best fit is obtained using a dust composition such that 95% of the grains are a mixture 
of amorphous silicates with the other 5% composed of iron oxide (17:1: 1 :1 ratio of interstellar 
silicates to glassy ol ivine to glassy pyroxene to wustite iDraine &: Led 1 1984 : iDorschner et al.l Il995l : 
Henning et al.lll995l ). We have extensively varied the grain comp osition using our lirn ited library 
and found this composition to yield the best result. Wustite (FeO: iHenning et al.lll995l ) is included 
as it provides a slight source of long wavelength emission between 40 and 50 jim which could not 
be fit by extending the outer shell radius or increasing the maximum grain size. However, even by 
including wustite, our model still has a slight deficit of emission in the 40 — 50 //m region. Wustite 
is featureless and the boost in emission is minor, thus we can not say with any certainty whether 
it is present in HD 179821. 

Our spherical shell result is presented in Fig. [H The model spectrum (solid dark line) is plotted 
on the Ay = 2 dereddened, T = 5750 K spectrum. Even though the central radiation source is a 
F-star (T = 7000 K), an optical depth through the shell of Ty = 1.95 obscures the photosphere 
so that it resembles a G-star (T = 5750 K). In particular the overall fit appears quite good in 
fitting data longward of 40 //m. Both our disc and shell models produce similar results for the 
ISO spectrum. Sub-millimeter and millimeter photometry highlights the main difference between a 
shell spectrum (small range of temperatures) and a disc spectrum (broad range of temperatures). 

The best shell fit yields an inner radius of = 9.8 x 10^ AU wh ich is in excellent agreement 
with mid-IR images and polarimetry o bservations (iGledhill et al.ll200lh . The placement of the inner 
shell radius is also in agreement with ISurendiranath et al.ll2002l . however we obtain an outer shell 
radius, ro, that is a factor of 6 smaller. This result may be more in line with polarization images 
in the near-IR which show a circul arly symmetric reflection neb ul a with a diaineter o f 15", a factor 
of ~ 0.75 times our outer radius (IKastner &: Weintraubl Il993l ). iBuemi et al.l 120071 assumed to a 
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n{a) oc dust distribution power law when fitting HD 179821 with DUSTY. While they arrived 
at a slightly cooler, closer shell, we found a satisfactory fit without appealing to steeper distributions 
by using a hotter central star. It is useful to m ention that we reproduced both previous model fits 
(jSurendiranath et al.ll2002l : iBuemi et al.l 120071 ). When the ISO spectrum was included, neither of 
those two previous models matched observations. In general, with a Tejj = 5750 K central source, 
we could not reproduce the infrared excess. To match the shape and peak of the infrared excess 
required a hotter, obscured central source (T^ff = 7000 K). 

Additionally from our fit, the dust temperature in the outer shell is T{ro) = 41 K. DUSTY also 
provides an estimate of the mass outflow rate although this value should be treated with caution as 
it has an inherent uncertainty of ~ 30%. Varying the gravitational correction by 50% produces no 
disc ernible effect on t he spectrum and is responsible for the uncertainty in the mass outflow rate 
(see llvezic et al.lll999l fo r more detail) . Our derive d outflow rate is M = 2 x 10 ~^ Mcr-, yr~ ^ 
in good agreement with ISurendiranath et al.ll2002l and a factor of 
We summarize results below. 



5 less then Buemi et al. 



agam 
2007i . 



6. Summary 

We have modeled the broad-band spectral energy distribution of HD 179821 using two distinct 
radiative transfer codes: one corresponding to a spherical shell geometry with p{r) oc r~^, and one 
to a disc geometry with p (r) oc r~^/^. Under these assumptions, both codes provided equally good 
fits for similar dust compositions and size distributions shortward of 40 //m. However, longward of 
40 /um, only the spherical shell model reproduces the observations. The radial dust density profile 
and corresponding range of temperatures present in the disc provide excess emission above the 
spherical shell solution and cannot fit the data. A wavelength of ~ 40 /xm marks the boundary for 
discerning the model fits. 

It should be noted that the density profiles need not uniquely indicate a particular geometry. 
If the densities deviate from the anticipated theoretical values of p{r) oc r^^ (shell), p{r) oc r^^/'^ 
(disc) , then the infrared excess is likely coupled to the density distribution rather than the geometry. 
The effect of changing the density profile in both codes should be investigated before comparative 
spectral modeling is established as a reliable method for distinguishing geometric differences. 

Based on our detailed spectral modeling, we conclude that it is likely that the nebula around 
HD 179821 is ci radicttivGly driven shell. Xlie dust is iiiciinly composed of sms^ll (^cirnax — 

0.25 pm) 

amorphous, glassy silicates. Interior to the shell is an evacuated region in addition to the central 
radiation source. The central star is most likely a T = 7000 K star obscured such that it appears 
as a T = 5750 K star (optical depth of ry ~ 1.95 through the spherical shell). This may aid in 
explaining previous conflicting spectral classifications. The dimensions of our shell agree well with 
previous observations and provide a good fit to the ISO spectrum. 

We compared the ISO spectrum to that of HD 161796, a confirmed post-AGB object (e.g. 
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Fig. [T]). In particular, the ISO spectra and broad-band SED look remarkably similar. Because our 
results scale with distance, our spectral modeling does not provide insight into whether HD 179821 
is a post-AGB or post-RSG. However, the similarity of the ISO spectra suggests that HD 161796 
and HD 179821 experienced a very similar mass-loss history. 
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Fig. 1.— Top: ISO SWS spectra of HD 179821 and post-AGB object HD 161796. In both objects, 
there is a steep rise between the 10 and 20 fim features possibly indicating a transition region to 
the opticahy thick outer wah or shell. HD 179821 was de-reddened using Ay = 2 while HD 161796 
was corrected using Ay = 1-2. Bottom: ISO spectrum of the post-RSG IRC +10420. 
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Table 1. Photometric Data 



Object 


U 


B 


V 


J 


H 


Ks 


HD 179821 


10.81 


9.49 


7.89 


5.371 


4.998 


4.728 



Note. — Magnitude summary for HD 179821. UBV mea- 
surei nents are on the Jo hnson system and are reproduced 
from^ Hrivnak et al. 1989 . JHKs measurements are fr om the 
2MASS point source catalogue (jSkrutskie et al.l 120061 ). 



Table 2. Colour Corrected IRAS fluxes, Submilhmeter Data and ISO Spectra for HD 179821 









FlOO/jm 


F45O Aim 




FlloO/^m 


ISO TDT 


32± 3 


700 ± 50 


430 ± 30 


160 ± 30 


1.22 ± 0.12 


0.227 ± 0.017 


0.058 ± 0.009 


11301444 



No te. — All presented fluxes have units of Janskys. Dat a is reproduced from Ivan der Veen et al 



19941 . ISO spectra is processed according to lSloan et al.ll2003l where multiplicative corrections are used 



for all wavelengths. 



Table 3. Model Summary 



Model 


Composition 


Teff 


To 




^max 




TV 5^ = 77 


Fits > 40 /im 


Disc 




K 


K 


/um 


^m 




AU 




Fig. El 


Pyroxene 


7000 


128 


0.05 


1.0 


0.25 


1.73x10^ 


No 


Fig. El 


Pyroxene 


7000 


128 


0.05 


1.0 


0.45 


1.69x10^ 


No 


Fig. El 


Olivine 


7000 


131 


0.005 


5.0 


0.25 


3.40x10^ 


No 


Fig. El 


Bronzite 


7000 


128 


0.005 


2.0 


0.25 


3.70x10^ 


No 


Shell 


Fig. H 


Mixture 


7000 


128 


0.005 


0.25 




1.95 15 


Yes 



Note. — 



Summary of our disc and spherical shell models. 
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A(/xm) 



Fig. 2. — Our wall models. The top figure corresponds to /x = 0.25 while the bottom corresponds 
to fi = 0.45. 
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A(/xm) 

Fig. 3. — Another fit to the wall model for glassy olivine (top) and glassy bronzite (bottom). The 
inclination in both figures is /x = 0.25. 
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X()um) 

Fig. 4. — Our spherical shell model (dark line). The dust temperature at the inner and outer radii 
are 128 K and 41 K. Even though the central radiation source is a T = 7000 K photosphere, the 
detached shell reddens the central source so that it resembles a T = 5750 K photosphere (thin line). 
The silicate features are fit using a mixture of glassy, amorphous silicates with a small component 
of FeO. 



